Effect of spin diffusion on spin torque in magnetic nanopillars 
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We present systematic magnetoelectronic measurements of magnetic nanopillars with different 
structures of polarizing magnetic layers. The magnetic reversal at small magnetic field, the onset 
of magnetic dynamics at larger field, and the magnetoresistance exhibit a significant dependence 
on the type of the polarizing layer. We performed detailed quantitative modeling showing that the 
differences can be explained by the effects of spin-dependent electron diffusion. 

PACS numbers: 72.25.Ba, 72.25.RB, 75.47.De 



According to the spin torque (ST) modey», current- 
induced magnetic switching (CIMS) in magnetic multi- 
layers is caused by angular momentum transfer from the 
conduction electrons to the magnetic layers. ST is be- 
lieved to occur within atomic distances from the mag- 
netic interfaces. Nevertheless, theories have shown that 
electron diffusion in the layers has an important effect 
on ST.— 1^1^ As a simple example, an electron scattering 
between two ferromagnets transfers angular momentum 
upon each reflection. However, this transfer is not neces- 
sarily associated with a net charge current /. Therefore, 
efficient utilization of electron scattering can result in re- 
duced I required to manipulate magnetic devices with 
ST. In a more subtle manifestation, spin-dependent elec- 
tron diffusion causes an asymmetry between the ST in 
antiparallel (AP) and parallel (P) configurations of the 
magnetic layers In an extreme case of such asymmetry, 
ST can change direction, resulting in anomalous current- 
induced behaviors?^ 

Despite extensive theoretical work, few experiments 
addressed the effects of diffusion on STi^ii^iaiia The 
main difficulty stems from the limited knowledge about 
the transport properties of individual layers in magnetic 
nanostructures. Different deposition and measurement 
techniques yielded significantly different values. On the 
other hand, both the Magnetoresistance (MR) and CIMS 
depend on the same spin-dependent transport properties. 
Therefore, simultaneous measurements of MR and CIMS, 
and their analysis within the same theoretical framework 
can lead to better understanding of the electron diffusion 
and its effect on ST. 

We report systematic measurements of MR and CIMS 
in nanopillar spin valves F1/N/F2 with identical free lay- 
ers F2=Py(5), Py=Ni8oFe2o, and different polarizers Fi 
incorporating Co. Thicknesses are given in nm. Large 
spin diffusion length Igf.Co makes Co ideal for studying 
the effects of diffusion. We used Fi=Co(20) in samples 
labeled Co20. To separate the contributions of the Co in- 
terfaces and its bulk, we tested samples with Fi=Co(3), 
labeled Co3, in which the scattering in the bulk of Co(3) 
was negligible. To eliminate the spin diffusion in the 
sample contacts, we inserted a strongly spin flipping bi- 
layer Fe5oMn5o(l)/Cu(l) between Co(3) and the bottom 
contact in samples labeled FeMnCo3. 

The multilayers Cu(50)/Fi/Cu(10)/F2/Cu(200) were 
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FIG. I: (a) dV/dl vs I at labeled Ji" and T = 295 K. Curves 
are offset for clarity, (b) same as (a), at T = 5 K. (c) J"*", I~ 
vs T for a Co20 sample, (d) Ic vs T measured ai H — 500 Oe 
for the three types of samples as labeled. 



deposited at room temperature 295 K (RT) by magnetron 
sputtering at base pressure of 5 x 10~^ Torr, in 5 mTorr 
of purified Ar. F2 and part of the Cu(lO) spacer were 
patterned into an elliptical nanopillar with approximate 
dimensions 130 x 60 nm. We measured dV/dl with four- 
probes and lock-in detection. Positive I flowed from Fi 
to F2. Magnetic field H was in the film plane and along 
the nanopillar easy axis. At least three nanopillars of 
each type were tested with similar results. 

Figs. [Ua),(b) show dV/dl vs I for a Co20 sample, ac- 
quired at RT and 5 K, respectively. The data at small 
= 50 Oe are characterized by hysteretic jumps to 
the P state with low resistance Rp at /~ < 0, and to 
the AP state with high resistance Rap at > 0. At 
H = 300/360 Oe in Figs. WLa)/{h), the jumps are re- 
placed by large peaks caused by the reversible transition 
between the P and AP states.^- The onset of the mag- 
netic dynamics starting at I = Ic appears as a sharp 
increase of dV/dl nearly independent of (1 kOe data 
in Figs. [Ua),(b)). The approximate equality Ic ~ I^ 
shown with a dashed line indicates that the reversal oc- 
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FIG. 2: (a) P-state resistances Rp offset by values at 5 K, 
and (b) MR va T for the three types of samples as labeled. 



curs when large-amplitude dynamics is excited by ST. 
The 5 K data exhibit significantly increased reversal cur- 
rents and Ic- Fig.[T{c) summarizes the temperature de- 
pendence of /+ and /~ . Both are nearly constant above 
130 K, below which they dramatically increase. Similar 
behaviors of Co/Cu/Co nanopillars indicate their intrin- 
sic origin from the spin-dependent transport in Co.^ 

One may attribute some of the dependence on T shown 
in Fig. [T]^c) to the effects of thermal activation. Indeed, 
J+ < Ic at RT because thermal fluctuations result in 
reversal slightly before the onset of large-amplitude dy- 
namics. In contrast, /+ > at 5 K because current- 
induced dynamics can occur before the reversal occurs. 
The fundamental quantity predicted by the models of ST 
is Ic- It is insensitive to thermal fluctuations and sam- 
ple shape imperfections, and can be directly determined 
from the sharp increase of dV/dl at H large enough to 
suppress hysteretic reversal. Fig. [Hd) summarizes Ic vs 
T for all three different sample structures. FeMnCoS 
data are approximately independent of T, while Ic for 
Co3 and Co20 increase when T is decreased. Compar- 
ing panels (c) and (d) reveals that Ic closely follows /+. 
It is not possible to measure a similar excitation onset 
current I^ in the AP state, because transition to the P 
state is not suppressed at any H. Below, we use /" as 
an approximation for I^ . 

Since F2 is identical in all samples, the different be- 
haviors of Ic must be attributed to the different spin- 
dependent transport properties of Fi. The difference be- 
tween Co3 and FeMnCoS is due to the spin flipping 
in FeMn, which eliminates spin diffusion in the bottom 
Cu(50) contact. The difference between the Co20 and 
Co3 data indicates that the effects of spin diffusion in 
Co are stronger than those in Cu. Despite a significant 
increase of Ic in Co20, it does not diverge as would be 
expected if the sign of ST was reversedi^ 

Figs. [Ha),(b) show temperature dependence of Rp — 
Rp{0) and MR= Rap — Rp- Rp increased with T due 
to magnon and phonon scattering, and were surprisingly 
consistent among the samples. Interestingly, there is a 
clear correlation between the variations of MR and Ic in 
all samples. As T decreases from RT, all MRs increase at 
a similar rate, while Ic slightly increase. At lower T, the 
trends for Co3 and FeMnCoS remain the same, while 



FIG. 3: (a) Resistivities of 40 nm thick Py, Co, and Cu films 
measured in Van der Pauw geometry. The Co and Cu data are 
fitted with the Bloch-Gruneisen approximation, with Debye 
temperatures 6co = 373 K and 6cu = 265 K. The Py data 
are fitted with a quadratic dependence, (b) Dependencies of 
residual resistivities on inverse film thickness (symbols) , with 
linear fits shown. 

a decrease of MR in Co20 at T < 130 coincides with a 
sharp increase of Ic- 

To understand the dependencies of MR, CIMS, and 
Ic on the sample structure, we performed simultaneous 
calculations of spin-dependent transport and ST. Our 
model combines a diffusive approximation for the ferro- 
magnets and outer sample contacts with a ballistic ap- 
proximation for the Cu(lO) spacer between the ferromag- 
nets,^ This approximation is consistent with calculations 
based on the Boltzmann equation.— We combine the con- 
tinuity conditions for spin currents and spin accumula- 
tion in the spacer between Fi and F2 derived by Slon- 
czewski^ (Eqs. (13), (14)) with a small-angle expansion of 
Eq. (28) for ST. The resulting expression for ST in terms 
of the spin current Ig — I^ — I^ and spin accumulation 
A/i = /i^ — /i-^ in the Cu(lO) spacer near the coUinear 
magnetic configuration is 

T ^ ^^iAGA^l - Is) (1) 

where e is the electron charge, Ti is the Planck's constant, 
G is twice the mixing conductance introduced in the cir- 
cuit theory,^ A is the area of the nanopillar, and 6 is 
the angle between the magnetic moments. At / = /c, r 
compensates the damping torque, yielding 

aejS22TTM2 , , 

= , (2) 

T 

where a ~ 0.03 is the Gilbert damping parameter^ 
7 is the gyromagnetic ratio, r is ST determined from 
Equation ([1]) at /= 1 in appropriate units, and S2 = 
M2V/2^B is the total spin of the Py(5) nanopillar. Here, 
V is the volume of F2, and /is is the Bohr magneton. 
The magnetization M2 of Py varied from 730 emu/cm^ 
at 20 K to 675 emu/cm^ at 300 K, as determined by mag- 
netometry of Py(5) films prepared under the same con- 
ditions as the nanopillars. These values are lower than 
expected for bulk Py, but consistent with the published 
results for Py films.— 
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FIG. 4: (a) Calculated Ic, Ic '"s G for FeMnCoS. (b) Same 
vs tcu for Co3. (c) Same vs lsf,Co for Co20, (d) same vs T 
for Co20 samples, for the residual values of lsf,Co as labeled. 



Equations H]) and ([2]) express Ic in terms of A/i and 
Is, the same quantities that determine MR in magnetic 
multilayers. We calculated Afi and Is self-consistently 
using a one-dimensional diffusive approximation employ- 
ing the standard MR parameters: spin asymmetries (3, 
renormalized resistivities p* = p/{l — 0^), spin diffusion 
lengths Igf in the layers, and similarly defined parame- 
ters AR* , 7, and 5 for the interfacesjii We estimate these 
parameters from a combination of the published values'-''' 
and our own measurements, as described below. 

The resistivity of each layer in our samples provides es- 
sential information about electron diffusion. Because of 
variations among published resistivities, we instead de- 
termined their values from measurements of thin films 
prepared under the same conditions as the nanopillars, 
with thicknesses verified by x-ray reflectometry. Fig. 
shows p{T) for 40 nm thick Py, Co, and Cu films, together 
with fittings for Co and Cu with the Bloch-Gruneisen 
approximation. We obtained better fitting for Py data 
with a quadratic dependence, indicating that electron- 
magnon scattering may dominate electron-phonon scat- 
teringji^ The dependence of the residual resistivity on 
film thickness was consistent with the Fuchs-Sommerfield 
approximation (Fig.[3l^b)), allowing us to extract the bulk 
residual values ppy(O) — 11.3 jUfJcm, pco(O) — 4.4 /xficm, 
and /Oc«(0) = 1.1 /iilcm. We used the extracted bulk 
p{T) to model all the extended layers in the nanopillars. 
The effect of lateral confinement in Py(5) nanopillars was 
approximated by using the resistivity of a Py(40) film. 

To estimate lsf{T), we used its empirical inverse re- 
lationship with p, along with the bulk residual values 
lsf,Py{Q) — 6 nm, and lsf,Cu(0) = 300 nm based on pub- 
lished measurements,— scaled by the somewhat different 
residual resistivities of our films. If scattering by thermal 
excitations does not flip electron spins, a weaker depen- 
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FIG. 5: (a) Calculated Ic vs T, and (b) calculated MR vs T 
for three sample types as labeled. 



dence lsf{T) oc ^Jl/ p{T) is possible. However, we show 
below that a dependence even stronger than l/p is more 
likely. We use (3py = ^iPy/cu = 0.7, ^ico/Cu = 0-8, /3co = 
0.36 for spin asymmetries, ^^Co/Cu ~ ^-^^ /^'^i^, 
^^Pj//c« ~ filrn^ for renormalized interface resis- 
tances, and Sco/Cu — 0-2, ^Pyjcu — 0-25 for spin flip- 
ping coefficients. ^■^'^ Their dependence on T is neglected 
due to the dominance of the band structure and impurity 
scattering far from the Curie temperature. For FeMn, we 
used lsf,FeMn ~ 0.5 nm, and ppeMn — 87 /ii7cm. Scat- 
tering at its interfaces was modeled by adding 0.5 nm 
to the nominal thickness of FeMn. To account for the 
Cu contacts, the calculation included outer Cu layers of 
thickness tcu, determined as described below. These lay- 
ers were terminated with fictitious spin sinks. 

To demonstrate that CIMS is extremely sensitive to 
the effects of diffusion, we now describe how our 5 K data 
can be fitted by appropriate choice of three parameters 
whose values have the largest uncertainty: conductance 
G in Equation ([T]), effective MR-active thickness tcu of 
the Cu contacts, and spin diffusion length lsf,Co- Calcu- 
lations for FeMnCo3 were significantly affected only by 
G, which controls the asymmetry of CIMS. The values 
of Ic/\I~\ in mA measured at 5 K for three FeMnCoS 
samples were 2.3/0.8, 1.6/0.6, and 3.1/1.5, giving an av- 
erage ratio Ic/\I~\ = 2.6. The calculated value increases 
from 1.46 at G = 0.5 fn-^m-^ to 6.1 at G = 2 fn-^m''^ 
(Fig.HIa)). The best values Icl\Ic\ = 3.34/1.27 are ob- 
tained at G = 0.87 frt~^'m~^, in reasonable agreement 
with band structure calculationsj^iiS, 

Spin diffusion in the bottom Cu layer has little effect 
on Go20 and FeMnCo due to the spin relaxation in Co 
and FeMn, respectively. To determine tcu, we use the ra- 
tios Ic/\I^\ of the three Go3 samples, 3.55/1.0, 4.6/1.5, 
and 4.2/1.2, giving an average ratio Ic/\I^ \ = 3.4. The 
calculated Ic/\Iq \ increases from 1.9 for tcu = to 14 
for tcu — 140 nm (Fig. IHb)), and eventually diverges at 
tcu = 200 nm. The best agreement with data is obtained 
for tcu = 55 nm, resulting in /c/|/^| = 4.4/1.3. 

Lastly, diffusion in Co significantly affects CIMS in 
samples Go20, but not in Go3 and FeMnCoi. We deter- 
mine Isf.Co irom the ratio Ic/\I~\ of five Go20 samples, 
8.9/2.1, 7.3/1.6, 9.0/2.0, 8.5/2.0, 8.0/1.7, giving an av- 
erage ratio Ic/\I~\ = 4.4. Fig. UJ^c) illustrates that the 
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calculated ratio Ic/\Iq\ increases from 1.0 for lsf,Co = 
to 5.2 for lsf,Co — 100 nm. The best agreement with the 
data is obtained for Igf.co — 42 nm consistent with the 
published values. 

Despite the ability to model the 5 K data, the cal- 
culations did not reproduce the dramatic dependence of 
Ic on T in Fig. [5] (see below). Therefore, one can at- 
tempt to determine lsf,Co from the dependence of Ic on 
T. Fig. [4l^d) shows calculations for the residual values 
hf,Co — 20 nm, 40 nm, and 60 nm. Large lsf,Co results 
in Iq decreasing with T, which is inconsistent with the 
data. Small lsf,Co gives decrease of Ic with T in better 
qualitative agreement with data, but gives unreasonably 
small /c at 5 K. Consequently, we return to the value 
determined from Fig. HJc) . 

Fig. EJa) shows the calculated Ic vs T for the three 
sample types. To interpret these results, we note that 
Figs. Iljb)-(d) exhibited an increase of Ic when the ef- 
fective MR-active resistance of Fi determined by pl^j 
was increased. This relationship was also established an- 
alytically^iSi The experimental correlation between the 
decreases of MR and increases of Ic in FigslU El is of 
the same origin. The lack of temperature dependence for 
FeMnCoS is therefore consistent with negligible spin dif- 
fusion effects in Fi. In calculations for Co3, the increase 
of Ic with T is caused by the increased contribution 
tculsf,Cu to the effective resistance of Fi. Calculations 
for Co20 show a competition between the contribution of 
the bulk Co resistivity, which increases with T, and the 
contributions from the Cu(50) layer and the outer Co/Cu 
interface, which decrease with T due to the increased spin 
flipping in Co. However, both Co3 and Co20 calculations 
do not reproduce the data, suggesting that the effects of 
thermal scattering should be re-examined. 

The calculated dependence of MR on T was in over- 



all agreement with data for Co3 and FeMnCo3, but did 
not reproduce the decrease at T < 130 K seen in Co20 
data (Fig. Eljb)). The calculations overestimated the val- 
ues, suggesting that our samples may be larger than their 
nominal size. However, this seems to contradict the cal- 
culated temperature dependence of i?p consistent with 
the data (not shown) , and the values of Ic that are larger 
than the measured 5 K values. This discrepancy can be 
reduced e.g. by decreasing lsf,Py, which results in a de- 
creased MR without significantly affecting CIMS. 

The failure of Co20 calculations to capture the de- 
crease of MR and the increase of Jc at T < 130 K in- 
dicates that lsf,Co decreases with T more rapidly than 
the accepted hf oc 1/p, resulting in the reduction of 
the effective MR-active resistance pcolsf.Co- One pos- 
sible mechanism for such a strong dependence may be 
electron-magnon scattering which can result in electron 
spin flipping without significant momentum scattering. 
We leave more detailed and perhaps alternative explana- 
tions to future studies. 

To summarize, we performed magnetoelectronic mea- 
surements of nanopillars with three different polarizing 
magnetic layers. The samples exhibited different current- 
induced behaviors, attributed to the spin diffusion in the 
polarizing layer. The calculations reproduced the lower 
temperature behaviors with reasonable values of trans- 
port parameters. However, temperature dependencies 
of magnetoresistance and current-induced switching indi- 
cate that the effects of thermal scattering on spin trans- 
port are more significant than presently believed. 
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